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Abstract 

The anomaly in the vacuum expectation value of the product of axial and two vector 
currents (AVV) in QCD is investigated. The goal is to determine from its value the 
vr*^ — > 27 decay width with high precision. The sum rule for AVV formfactor is studied. 
The difference f^o — /^+ caused by strong interaction is calculated and appears to 
be small. The — t] mixing is accounted. The vr*^ — > 27 decay width determined 
theoretically from the axial anomaly is found to be r(7r'^ — > 27) = 7.93eV with an error 
~ 1.5%. The measurement of the 27 decay width at the same level of accuracy 

would allow one to achieve a high precision test of QCD. 

PACS numbers: 11.15.-q, ll.30.Qc, 12.38.Aw 
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The statement that in massless QED the axial anomaly is contributed by massless pseu- 
doscalar state had been first done by Dolgov and Zakharov in 1971 pj. Later this problem 
was investigated in many papers. (See, e.g. the book [2] and the reviews [3], [1]). Now it 
is known, that the anomaly for transition of axial isovector current into two massless vec- 
tor currents, i.e., into two photons, is almost completely exhausted by the contribution of 
7r° — > 27 decay. However, the accuracy of the calculation of the 7r° — 27 decay contribution 
is about 5-7% [3] and of the same order is the experimental error in the determination of 
the 7r° 27 decay width r(7r'^ —>■ 27) [5j. The PriMex experiment at JLab [6], where it is 
planned to reduce the error in the experimental value of r(7r° 27) down to 1-2%, is under 
way. Also, it is desirable to have theoretical prediction at the same level of accuracy and 
therefore, to have a possibility of high precision test of the anomaly -a significant ingredient 
of the modern field theory. 

The calculation of hadronic contribution to the AVV anomaly at the desired level of ac- 
curacy was done by Moussalam [7] and by Goity, Bernstein and Holstein [8] in the framework 
of the Chiral Effective Theory (GET) (in the second order) and of 1/Nc expansion. In these 
calculations the result is expressed through the parameters of the GET Lagrangian in the sec- 
ond order, which were determined from the set of the data. In the present paper we perform 
such calculation in QGD using the dispersion relation representation for the AVV formfactor, 
QGD sum rules to determine /^o — /^+ difference and accounting for tt^ — r] mixing. The only 
parameter, which enters the result, is the value r{ri ^ 27) 

The notation 

Tf,at3{p,p') = {P,ea;P,e'f^\jj^^\0) (1) 

is used for the matrix element of the transition of the isovector axial current 

jjS = {uit,l5U - (J7m75C?)/V2 (2) 

into two photons with momenta p, p' an polarizations ea , e^. (Here u and d are the fields of 
u and d quarks.) The general form of T^ai3{p,p') is (g = p + p') [9].[T0].p]: 

Tf^a|3{p,p') = Fi{q^)q^eal3paPpP'a + 

+ ^^2(g^)(Pae^/3pa - P'p^papa)PpP'a (3) 

The functions Fi{q^) and F2{q^) can be represented by nonsubtracted dispersion relations 
and in QGD the anomaly condition can be written as the sum rule [4j: 

oo 

j ImF,{q')dq' = V2a{el-el)N,, (4) 

where Cu and ^ are w-and d-quarks electric charges, e„ = 2/3, = —1/3, Nc is the number of 
colours, Nc = 3. Note, that at large the function ImFi{q^) ~ {l/q^)lnq'^ and the integral 
in (jl]) is well converging. Emphasize, that in QGD there are no perturbative corrections to 
Eq-dlJ [H] and it is expected, that nonperturbative corrections are absent also. 

Let us saturate the left-hand side (l.h.s) of (jlj) by the vr" contribution. Use the relation 

(0|jSV°) = (5) 
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The general form of the vr^-contribution to T^ai3{p,p') is 

Tf,af:i{P,P') = -/ttO^-^ — ^A^qf,eai3XaPxp'a (6) 
q 

where is a constant. In the approximation, when only pion contribution is accounted in 
the l.h.s. of (jl]) the constant is found by substituting into (jl]). The result is (numerical 
values of Cu, e^, Nc were used): 

AS> = - f . (7) 

The 7r° — > 27 decay width is easily calculated from ([6]), ([7]), 

Index (0) at F means that the saturation of the sum rule (jl]) by the vr^ state was exploited. 
At /^o = = 130.7 ± 0.4 MeV and m^o = 135.0 MeV |5] we get from ([8]): 

r(o)(7r" ^ 27) = 7.73el^ (9) 

Turn now to the calculation of correction to zero approximation. The first correction is 
due to the fact, that generally /^o is not equal to f^+. There are two sources of the /^o — 
difference: the electromagnetic interaction and violation of isospin in strong interaction, 
vr^ has no electromagnetic interaction and the electromagnetic interaction of Tr^ was already 
accounted, when the value of f.,^+ was found from the vr"*" fi^u decay data [5j. So, the 
electromagnetic interaction does not change the value of 7,^0 in comparison with the pre- 
sented above value of f^. (The discussion of magnitude of f^^o — difference, caused by 
electromagnetic interaction was done in Ref. [8]. In order to find /t^o — caused by strong 
interaction consider the polarization operators n|fj and U^^^ of axial currents jj^^ ([2]) and 

i^s^ = ^Itil^d. (The presented below method is exposed in [12], [13] and reviewed in [H]). 
The general form of these polarization operators is 

n(:)(g) = n«(g2)(g,g. - 5,,q') + q,q.Ilf{q'), ^ = 3, + (10) 

The pion contribution to n^,^(g) is given by [T2] : 

^^lu{.q)^ = — liQf^qu - 5^uq^) fq^q^, ^ . (11) 

Consider the longitudinal polarization operator 11^'' (g^) to which only pseudoscalar mesons 
are contributing. In order to separate the interesting for us second term in ffTTl) let us multiply 
ffTTj) by q^qu/q'^ and consider the difference 



nf (g^) - nW(g^)) = + (12) 

J q^ - m\ q^ - mU 



As was demonstrated by Weinberg [in], in the first order in 171^ — ma the 7r+ and 7r° masses are 
equal and the experimentally observed vr"*" and vr" mass difference arises from electromagnetic 
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interaction. In the second order in u- and d-quark masses m^o — fn'^+ is proportional to 
(m„ — nid)'^ or (m^ — rud) [(OImmIO) — (0|Jc?|0)] and, as can be shown, is very small. So, in 
(fT2|) we can put = m^o and have 

A(g^) ^ g^(nf (g^) - nW(g^)) = -^^{f.. - /^). (13) 

Exploit the standard QCD sum rule technique and represent A(g^) as an operator product 
expansion (OPE) 

A(g2) = i?2(g') + R,{q') + i?6(g'), (14) 

where R2 corresponds to the contribution of bare loop diagram, i?4 - to the operator of 
dimension 4 and Rq - to the operator of dimension 6. Only the terms, proportional to 
m„ — TUd and (OImmIO) — (0|(i(i|0) remain in f fT4l) . Note, that A(g^) is even under interchange 
u ^ d. For this reason no linear in m„ — md or (OImmIO) — (0|(i(i|0) terms can survive in the 
right-hand side of (IT^ . Calculating the terms of OPE we get after Borel transformation: 

md — rriu 

[au-ad), (15) 



where 



47r2M2 



-(27r)2(0|gg|0), q = u,d (16) 



and is the Borel parameter. The contribution of the d = 6 operator vanishes. For the 
difference /^o — we have 



A/. _ Uo - U+ 3 M2 



(m„ - md) - - — — — 70, 



(17) 



where 



1= (1^ 

ttq 



For numerical estimations we put — 771^ = 3.5MeV ± 20% [H], 7 = 6 ■ 10 ^ [16] or 
7 = 3- 10~^ fTT] and M = IGeV. The calculations give 

^-2-10-^ dS]; ^^4-10-^ [H] (19) 

In In 

Therefore, the difference /^o — /^+ is negligible. 

The additional source of the /^o — /^+ difference arises because there is the contribution 
of //-meson to the polarization operator besides the rr^. This contribution appears from 
the r] — 7r mixing and must be subtracted from the value of /^o. Omitting the details of simple 
calculations, we present the result 



(20) 
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-1.3.10"'' (21) 



where Orjir is the t] — it mixing angle given below (Eg. (1291) ). Numerically 

m 

and is also very small. 

Let us dwell now on the main correction to the zero approximation result for r(7r'' 27) 
arising from excited states contributions to the sum rule (jl]) . The next after vr" pseudoscalar 
meson is 7]. It can contribute to (jl]) because of isospin violation, caused by different masses 
of u- and d- quarks and resulting in 7r° — 77 mixing. The problem of 7r° — r/ mixing in QCD 
was considered in [18], [19], [20]. Following [18], [20] introduce nonorthogonal states IP3) and 
IPs) and the corresponding fields ^93, coupled to jj^^ and jj^^: 

{0\jlS\Pi)=^Skifkqf., k = 3,8. (22) 

Nonorthogonality of the fields ips, ips corresponds to the non-diagonal term AH = rn^^ipsfs 
in the effective interaction Hamiltonian. In the presence of such term the standard PCAC 
relations are modified in the following way [18\ 



dt^jjS = fr,{ml(f8 + ml^fs)- (23) 
The fields (^3, ips are expressed through the physical fields (p.„, as 

ips = -Sin9r,nV7T + Cos9^^(pr„ (24) 
where the mixing angle is given by 

0,. - - (25) 

In QCD the nondiagonal mass m^^ is expressed through the difference of u- and d-quark 
masses [T8].[T9]: 

= ^ 

V3 m„ + nid 

Now Im Fi{q^) is given by the sum of contributions of vr^ and 77 mesons. In order to separate 
the formfactor Fi{q^) and to kill the contributions of axial mesons in the representation of 
Im Tfj,ai3{p,p') in terms of physical states, multiply Eq.([3]) by q^/q"^- Using Eq.'s f|T]). fl23l) .f l24|) 
and taking the imaginary part, we get: 

1 1 

Im 9^^(27 I jjS I 0) = -^^"^(27 I mKCosOrinVn + Sinerjn^r,) + 

+ml^{-Sin9r,jr'-Pn + CosOnnipyj) I 0) ^ 7r[ 6{q^ - ml)A^+ 
■ <- 5(g^ - ml) A, - ^^^5{q' - ml)A^ + ^5(g^ - m^) A, ] ^ 



1 - 



m; 1 / niu - rrid 



+ 



1 



rrid 



mu + TUd 



\5{q^ - ml) A. 



V ( ' 



(27) 



where A,^ is the amphtude 77 — > 27 decay. The ratio Arj/A^, is equal 



A^ 

AjY 



r(?7 27) 
r(7rO ^ 27)^ 



1/2 



(2^ 



Numerically, at T{ri 27) = 510e1/ [5] and 1(71° ^ 27) = 7.73eV, A^^/A^ = 1.0. According 
to fl25l) . fl26p the mixing angle ^^tt is given by (the numerical data were taken from Ref.[14j): 



1 m„ 



rj-K 



y/Smu + rrid ml 



-0.0150 ±0.020 



(29) 



and the total correction to the amplitude of vr*^ 



27 decay arising from 77 — 71 mixing is 
-(1.2 ±0.25)%. (30) 



Let us estimate the contributions to the sum rule (jlj) from higher pseudoscalar mesons. 
The rj' meson contributes through rjri' mixing. Its contribution is suppressed by the ratio 
{jriri/mrfY ill comparison with rj contribution and we have 



nv' ^ 27) / 



3/2 



(Fi), \ T{rj ^ 27) \m,^ 



mr. 



m^ji 



Sind 



rjr]' 1 



where t?^^' is the 7777' mixing angle. In accord with ^21j we can put ^rfq' 
r(?7' — > 27) = 4.29 MeV [5J we have the estimation 

(Fi),V(Fi), ~ 0.15 



(31) 



-20°. Using 



(32) 



Finally, estimate the contribution of the resonance 7r(1300). The contribution of this reso- 
nance to Tf^af3{p, p') ^ (as well as any other excited states) should be proportional to mu+md 
or, what is equivalent, to m^, otherwise the axial current would not conserve in the limit 
m„ ± md 0. Therefore, for dimensional grounds. 



(1300) 



mz 



(Fi)^ 



< 1%. 



m 



(33) 



7r(1300) 



In fact, probably, fl33|) overestimates the 7r(1300) contribution, since, as was mentioned above, 
Fi(g^) decreases as l/q"^ at large q^. So, in what follows we take (-F'i)7r(i3oo)/(-Fi)7r ~ 0.5%. 
Summing all uncertainties in squares we have the total uncertainty in the amplitude about 
0.7%. The final result for r(7r° 27) from the axial anomaly is (the correction, given by 
Eq. (l2T!) is accounted): 

r(7r° ^ 27) = 7.93eV ± 1.5%. (34) 

In the limit of errors this agrees with those found in and [8]. The experimental test of 
this relation, which is planned by PriMex Collaboration would be very important: it would 
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be a high precision test of QCD and even more general - the phenomenon of the anomaly. 
Not too many such high precision tests of QCD are known till now. 
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